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ABSTRACT 



We extend our recent work, in which the Dirac equation with a "(asymp- 
totically free) Coulomb + (Lorentz scalar 7ocr) linear " potential is used to 
obtain the light quark wavefunction for qQ mesons in the limit mq — > oo, 
to estimate the decay constant fp and the mixing parameter B of the pseu- 
doscalar mesons. We compare our results for the evolution of fp and B with 
the meson mass Mp to the non-relativistic formulas for these quantities and 
show that there is a significant correction in the subasymptotic region. For 
a = 0.14 GeV -2 and A^ = 0.240 GeV we obtain: f D = 0.371 , f Ds = 
0.442 , f B = 0.301 , f Bs = 0.368 GeV and B D = 0.88 , B Ds = 0.89 , B B = 
0.95 , B Bs = 0.96 , and B K = 0.60. 



Recently, we presented our results for the form factor (Isgur-Wise function) which 
parametrizes the transition matrix elements of mesons containing a heavy quark [|IJ . In 
our relativistic treatment of a heavy-light system, we assumed that, with the heavy quark 
fixed at the origin, the light quark wavefunction obeys a Dirac equation with a spherically 
symmetric potential. This potential, an asymptotically free Coulomb term plus a linear 
confining term, reflects the QCD interactions both at short and long distances. In 
this way, we obtained the shape and the slope at zero recoil of the form factor in the 
leading 1/w,q order by fitting our model parameters to the experimental results on the 
semileptonic B decays. 

In this paper, we extend our results by calculating the decay constants fp and mixing 
parameters B of the pseudoscalar qQ mesons in our model, fp parametrizes the matrix 
element for the decay of a pseudoscalar meson through the axial- vector current A^: 

<0\Al(x)\P(k)>=tk,f P e- lk - x , (1) 

where k is the four-momentum of the meson. On the other hand, B is related to the 
matrix element for neutral-meson mixing which is conventionally written as 

< P|(V M - 4)> >=ppM P B , (2) 

where is the vector current and Mp is the meson mass. From a phenomenological 
point of view, knowledge of these matrix elements is necessary for extracting important 
quantities such as the quark mixing matrix and CP violation within the Standard Model. 
Experimentally f n = 132MeV and fx = 167MeV are well-known. However, for heavy 
mesons only a few data with large uncertainties are available. Mark III sets an upper 
bound f D < 290MeV § while WA750, CLEOIIgj and BES§ find 225 ± 45 ± 20 ± 
41 , 344 ± 37 ± 52 and 430±i|g ± 40MeV for f Ds respectively. 

There are various quark-model, QCD sum rules and lattice calculations of these 
parameters]!], [7|, ]9|]. A frequently used point of comparison is a scaling law for fp that 
is derived from a non-relativistic (NR) quark model: 



/^V Q -oo)ay^, (3) 

where tuq is the mass of the heavy quark. The mixing parameter is identically 1 in any 
NR quark model: 

B NR (m Q ) = 1 . (4) 

In Refs. H [J], it is shown that the combination of a relativistic dynamics and an 
asymptotically free Coulomb interaction for qQ system results in a significant deviation 
from the NR scaling law. Our approach here is complementary to Refs. [§, || in the 
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sense that we use a potential which is not only asymptotically free but also includes a 
linear confining term which determines the global behavior of the Dirac wavefunction. 
At the same time, we use a saturation value for the strong coupling a 
therefore avoiding unphysical pair creation effects. 
We start with the time-independent Dirac equation: 



oo 
s 



a.Jr 



oo 



-*V) + Vc(r) + c + 7°(ot + m o)\ = e q^(r) 



(5) 



where m q is the light quark mass. V c (r) is the asymptotically free Coulomb potential 
(tranforming as the zeroth component of a Lorentz four- vector) : 



V c (r) = -- 



(6) 



where a s (r) is obtained in the leading log approximation and is parametrized as follows 



2tt 



a.ir) 



;ii-M*)H5 + f] 



(7) 



The parameter 6 defines the "long distance" saturation value for a s . We use 5 = 2.0 
which corresponds to a s (r = oo) ~ 1.0. The parameter (3 is related to the QCD scale 
A-jjg by (3 = (2.23A]^)~ for Np = 3 (we use Np = 3 throughout this paper). We take 
Ajfg ~ 0.240 GeV , which corresponds to (3 = 1.87 GeV -1 in most of the calculations. 

To describe the long distance behaviour, a linear term (7oot) is introduced in the 
potential which transforms as a Lorentz scalar. This form is favored based on many 
theoretical and phenomenological arguments [K], 11, |I2 
mainly use a value a = 0.14 which is favored by a recent lattice estimate 



For the parameter a, we 
and can 
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also be extracted from the Regge slope using two-body generalization of Klein-Gordon 
equation [jlT], [T2| . However, we also present results for a = 0.25 and 0.18 GeV 2 , as the 
former value is compatible with the experimentally available 0.45 GeV 2P-1S splitting 
for D and D s system and the latter is extracted from Regge slope data using string 
model [13]. As mentioned in the discussion of our results in Ref. []T|, the best fit of our 
model to the recent CLEO 1994 data analysis favors a smaller value for the parameter 
a. 

The potential in the Dirac equation includes an additive constant cq, which is clearly 
subleading both in the short and long distance regimes. The only role of this constant 
is to define the absolute scale of the light quark energy e q which is identified with the 
"inertia" parameter often introduced in heavy quark effective theory: 



lim (M P - m Q ) 
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Therefore, only the difference A q — cq ~ e q — cq can be extracted from Eq. (5). As we 
indicate later, in the heavy limit mq — > oo, fp and B are not very sensitive to A q (and 
in turn to Co). 

The ground state solution to Eq.(5) has the form: 



*<f) = N ?'Lf-\ > ( 9 ) 




where iV is a normalization constant. Figure 1 illustrates the functional behaviour 
of the normalized large component x n (= Nx{r)) and small component n (= N<f)(r)) 
of the wavefunction ip(r) for the cases where m q = (appropriate for q = u,d) and 
m q = 0.175 GeV (for q = s). We observe that x{ r ) is finite and 0(r) — > as r — > 0. 
This, as noted in Refs. || is due to using a running a s (r) in the Coulomb poten- 
tial (Eq.(6)) rather than a fixed a which would result a singular wavefunction at the 
origin. On the other hand, since our model incorporates the long distance behaviour 
of QCD interactions (through linear confining potential) and at the same time avoids a 
singularity in the Coulomb potential (by introducing a saturation value for a s (r)), the 
resulting wavefunction is physically meaningful for the whole range of r. This is our 
main improvement over previous works where a Dirac equation along with the leading- 
log Coulomb potential has been applied to heavy-light mesons || [|. The functions x( r )> 
<p(r) and the normalization constant iV are independent of tuq in our leading order (in 
I/tuq) approach where we assume that the heavy quark is fixed at the origin. How- 
ever, as it is illustrated in Figure 1, the SU(3)f symmetry breaking strange quark mass 
m s = 0.175 GeV results in about 20% larger value for iV (x{r) is independent of m q as 
r — > 0). The normalization constant iV depends on the global behaviour of the poten- 
tial, i.e. P and a (see Eqs.(5), (6) and (7)). We mainly use (3 = 1.87 (corresponding to 
Ajjs = 0.240 GeV) and a = 0.14 GeV -1 (see the discussion following Eq. (7)), however, 
we do vary these parameters to investigate the sensitivity of our results. 

We now proceed to compute the decay constant fp and mixing parameter B in our 



model. The decay constant is given by the overlap integral |T5 

fp(m Q ) = ^~ p N 2 J d?r\^ Q {r)\\\r) , (10) 

where 12 is a color-spin-flavour coefficient and the quantum-mechanical fluctuations 
in the position of the heavy quark are assumed to be described by a non-relativistic 
wavefunction tyqif). Considering that ^q{t) is "spread" over a distance rq = 0(l/m<g), 
one can reduce (10) to the expression f9j: 



fp(m Q ) = J-^N X (r Q ). (11) 
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At this point, we need to make a definite connection between tq and the heavy quark 
mass rriQ. It is assumed that rq = k/ttiq, where k > 1 is expected on physical grounds. 
We obtain a value for k by extrapolating the decay constant formula (11) to K-meson 
system where fx = 167MeV is known experimentally This extrapolation, even though 
of uncertain reliability due to finite rriQ effects, has been frequently made in purely NR 
models. For constituent ~ m k = 0.495 GeV, we obtain k = 1.67 for = 0.240 and 
k = 1.73 for Aj^ = 0.200. In comparing our results for the decay constant with the 
NR scaling law, Eq. (3), we also include a finite renormalization factor suggested by 
Voloshin and Shifman [Jl|J and Politzer and Wise [1? 
between the "improved" NR scaling law 



In Figure 2, comparison is made 
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and our relativistic results 
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(12) 



(13) 



where in the latter we have made the approximation Mp « itlq in the expression for rq. 
One can justify this approximation in view of the uncertainty in the determination of 
tq from the extrapolation to K-meson system. Because of this assumption, rq does not 
depend on the inertia parameter A q and therefore is also independent of the constant 
c . Also in Eq. (13), 7 = 2/9 for N F = 3 and a s (M P ) = 2ir/9\n(M P /A m ) in the 
renormalization correction factor. In the Figure, the above functions are normalized 



to their value at a large mass scale 



60A 



MS* 



where the physics (of rriQ — > 00) 



is 



well understood. Our results are presented for m q = and m q = 0.175 GeV. The 
main feature that distinguishes the relativistic graph from the NR scaling law is the 
maximum at around D meson mass, as noted in Refs. |8], |9|. We would like to emphasize 
again, that our improved model is free of unphysical behaviour for the whole range 
of r. Therefore, here we rigorously confirm that the maximum in the scaling law is 
physical and is due entirely to the relativistic dynamics of the light quark at the short 
distances. The position of this maximum is around 7.0Aj^- for m q = (for = 
0.240 GeV , M£ eak w 1.68 GeV which is somewhat below M D w 1.87 GeV) and around 
8.lAjjg for m q = 0.175 GeV (resulting M|f ak « 1.94 GeV which is roughly the same as 
Md s ~ 1-96 GeV). From Figure 2, we also notice that, roughly speaking, from the B 
meson mass scale onward, SU(3)p is a good symmetry as far as the scaling of the decay 
constant is concerned. 

The mixing parameter B can also be written as an overlap integral 



B(m Q ) 



12 



fpMf 



N 2 d 3 r\V Q (r)t 



(14) 
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which will be approximated by 

(15) 

by the same argument that followed Eq.(lO). The non-relativistic scaling law for mixing 
parameter (Eq. (4)) is recovered for tjiq — > oo (i.e. tq — > 0) as 0(r — ► 0) = (see Fig. 

!)• 

Figure 3 illustrates the evolution of the mixing parameter B with the meson mass 
Mp (as in the case of the decay constant, tq = k/ttiq and we use the approximation 
Mp « ttiq) for m q = and m q = 0.175 GeV. We observe that in the subasymptotic 
region, deviation from the NR scaling law is 5% and 12% for the B and D meson mass 
scale respectively. On the other hand, the SU(3)p symmetry breaking effects in the 
scaling;, even for M D ~ M P ~ M B , is negligible. 

In Table 1, we present our numerical predictions for the decay constant and mixing 
parameter of D, D s , B, B s mesons plus an estimate of Bk, the mixing parameter for the 
K-meson. To test the sensitivity of these predictions, we vary the model parameters (3 
and a. For example, changing from 0.240 GeV {(3 = 1.87 GeV" 1 ) to 0.200 GeV 
((3 = 2.24 GeV -1 ), results in a very small change (ranging from 1% to 4%) in fp and 
B. The ratios /b s //b = 1-22 and foj ' fo = 1-18 (for a = 0.14 GeV -1 ) are independent 
of A^r. On the other hand, the sensitivity to the linear potential coefficient is much 
more significant. Besides a = 0.14 GeV -1 , the decay constants and mixing parameters 
are also estimated for a = 0.18, 0.25 GeV -1 (see the paragraph on the linear potential 
following Eq. (7)). An increase of 30% to 40% in the decay constants is observed for 
increasing a from 0.14 to 0.25 GeV -1 . However, the mixing parameters are far less 
sensitive to this model parameter. 

The results obtained for f D all exceed the upper bound of 0.290 of Ref. 2. However, 
the finite renormalization factor is uncertain at the lower mass of the D meson; if this 
factor, which is 1.27, is not included, fp, is 0.292 for o = 0.14. The results may therefore 
not be incompatible with Ref. 2 at the smaller value of a. 

A comparison between our predictions (Table 1) and other theoretical results || also 
favors a smaller value for o. Even though our results are larger than other theoretical 
predictions, for a = 0.14 GeV -1 and taking into account the uncertainty factor 1.27 
(from the finite renormalization factor below D-meson mass scale) there are agreements 
with some lattice and potential model estimates. However, QCD sum rules predictions 
are generally smaller than ours. Our estimated Bp and Bp s is on the low side of a recent 
lattice prediction by UKQCD0. 

In conclusion, we used a relativistic model with a phenomenological potential that 
accounts for QCD interactions at all length scales, to estimate the decay constant fp 
and mixing parameter B of heavy-light mesons. The evolution of fp and B with the 



B(m Q ) = 1 - 
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meson mass significantly deviates from the NR scaling law in the phenomenologically 
interesting subasymptotic region M D ~ M P ~ M B . 
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Table Captions 

Table 1: The estimated decay constants (in GeV) and mixing parameters for 
various choices of the model parameters. 
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Figure Captions 



Figure 1 '. The normalized large \ n = ^x{ r ) an d small <p n = N<p(r) component of 
the wavefunction for q = u,d (m q = 0) and q = s (m q = 0.175 GeV). 



Figure 2: The evolution of the decay constant with the meson mass (m q = and 
m q = 0.175 GeV) compared with the NR scaling law. 



Figure 3: The scaling of the mixing parameter B with the meson mass. 
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Figure 3 
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Table 1 





a = 0.25 


a = 0.18 


u = 0.14 


P = 1.87 


f D = 0.519 f B = 0.427 
f Ds = 0.587 f Bs = 0.496 
B D = 0.87 B B = 0.95 
5 Ds = 0.89 B Bs = 0.95 
B K = 0.58 


/d = 0.429 f B = 0.350 
/d, = 0.497 f Bs = 0.417 
B D = 0.87 5 B = 0.95 
B Ds = 0.89 5 Bs = 0.95 
B K = 0.59 


f D = 0.371 / B = 0.301 
f Ds = 0.442 / B< = 0.368 
B D = 0.88 B B = 0.95 
B Ds = 0.89 B Bs = 0.96 
5^ = 0.60 


(3 = 2.24 


/d = 0.514 f B = 0.418 
/d. = 0.580 f Bs = 0.485 
B D = 0.88 5^ = 0.95 
B Ds = 0.89 5 Bs = 0.96 
B K = 0.59 


f D = 0.424 / B = 0.341 
f Ds = 0.491 / Bs = 0.406 
B D = 0.88 B B = 0.95 
5 Ds = 0.90 B Bs = 0.96 
5 X = 0.60 


f D = 0.367 f B = 0.292 
/ Ds = 0.434 f Bs = 0.356 
Bb = 0.88 B B = 0.95 
B Ds = 0.90 5 Bs = 0.96 
B K = 0.61 
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